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. In this paper, we examine in a unified fashion dissipative transport in strongly correlated systems. We thereby 

demonstrate the connection between "bad metals" (such as the high temperature superconductors) and "perfect 
fluids" (such as the ultracold Fermi gases, near unitarity). One motivation of this work is to communicate to the 
high energy physics community some of the central unsolved problems in high T c superconductors. Because 
of interest in the nearly perfect fluidity of the cold gases and because of new tools such as the AdS/CFT corre- 
spondence, this better communication may lead to important progress in a variety of different fields. A second 
motivation is to draw attention to the great power of transport measurements which more directly reflect the 
(— ( , excitation spectrum than, say, thermodynamics and thus strongly constrain microscopic theories of correlated 

q ■ fermionic superfluids. Our calculations show that bad metal and perfect fluid behavior is associated with the 

' presence of a normal state excitation gap which suppresses the effective number of carriers leading to anoma- 

lously low conductivity and viscosity above the transition temperature T c . Below T c we demonstrate that the 
condensate collective modes ("phonons") do not couple to transverse probes such as the shear viscosity. As a 
result, our calculated shear viscosity at low T becomes arbitrarily small as observed in experiments. In both 
ryj , homogeneous and trap calculations we do not find the upturn in r\ or r//s (where s is the entropy density) found 

. in most theories. In the process of these studies we demonstrate compatibility with the transverse sum rule and 

' find reasonable agreement with both viscosity and cuprate conductivity experiments. 

&' 

PACS numbers: 03.75. Ss,67.10.Jn, 67.85.De 
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I. INTRODUCTION 

l o j . 

There is a widespread interest in studying ultracold Fermi gases to learn about some of the most strongly interacting systems in 
nature. It has been argued 1 that the atomic Fermi gas near unitarity provides a prototype for the physics of quark-gluon plasmas 
which are associated with the early stages of the Big Bang. Most remarkable about these two nominally different physical 
systems is the fact that they correspond to nearly perfect fluids, exhibiting even in their non-superfluid phases extremely small 
values of the shear viscosity rj. In this paper we call attention to a third class of strongly interacting systems which is the metallic 
counterpart of the perfect (neutral) fluid, namely, the so-called "bad metal" 2 . We show how low conductivity and low viscosity 
, are analogous. Bad metals are known to exist in nature in a most fascinating class of materials: the high T c superconductors. 
We argue here that to learn more about the physics of nearly perfect fluids it is particularly useful to study the behavior of the 
conductivity in the cuprate superconductors and vice versa. We note that it may seem at first sight paradoxical that bad metals 
t— I . and perfect fluids have anything in common since the condensate contribution to the conductivity is infinite in the metallic 
system whereas the viscosity counterpart in the neutral superfluid is zero. We stress throughout this paper that the dc transport 
• i-h , (associated with ui — > 0) in both cases reflects the normal fluid or excitations of the condensate. 

Because the ultracold Fermi gases near the unitary limit are also thought to be related to quark-gluon plasmas 1 , much attention 
has focused on particle-physics-based calculations of the anomalously low shear viscosity, ?f' 4 . However, considerable insight 
into the thermodynamics 5 and various spectroscopic studies 6 7 has also been obtained via a condensed matter perspective. This 
paper belongs to the second school in which BCS theory is extended to accommodate arbitrarily strong interactions. We apply 
BCS to Bose Einstein condensation (BEC) crossover theory 6 to compute r\ in the neutral Fermi gases and the conductivity a 
(evaluated at to — > 0) for the charged counterparts, such as the high T c superconductors. Our results, which are reasonably 
consistent with experiments, apply both above and below T c . Essential to this work is demonstrable consistency with central 
sum rules and conservation laws 8 . 

Experimentally this crossover can be straightforwardly observed in the ultracold Fermi gases by exploiting Feshbach reso- 
nances. As has been argued 9 , this crossover also appears relevant to the cuprates (i) because of their anomalously short coherence 
length. Added support comes from (ii) their anomalously high transition temperatures and (iii) the widespread effects of a normal 
state gap 6 ' 10 , (iv) The smooth evolution from over to underdoped behavior of the optical conductivity, has similarly been used to 
argue 1 1 that the BCS-BEC crossover scenario might be appropriate to the cuprates. Although the cold gases and high T c cuprates 
have s and d-wave order parameters respectively, this difference does not significantly modify our treatment of transport. 
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The shear viscosity, like the dc conductivity, is a powerful probe for testing microscopic theories because it reflects the 
normal fluid component and is extremely sensitive to the nature of the excitation spectrum. The low T normal Fermi liquid 
phase with scattering lifetime 7 _1 and effective mass m* helps to set up the general nomenclature. Here r\ — |m;fi7 _1 m*. 

2 _ , 

Similarly a = ^fr7 ■ More generally, one can think of r\ and a as characterized by the effective number of normal excitations 
(n — > n e f f (T)) as well as their lifetime which we emphasize here is a many body effect. Crucial is an understanding of how 
n e f f depends on T. 

Our central conclusion is that both the effects of a fermionic gap (with onset temperature T* > T c ) and non-condensed pairs 
act in concert to reduce n e // and thus lower the shear viscosity and dc conductivity at all T < T*. These non-condensed 
pairs are associated with the stronger than BCS attraction and are present for < T < T*. When quantitatively compared 
with very recent shear viscosity experiments 12 (we independently infer 13 an estimated lifetime from radio frequency data) the 
agreement is reasonable. We similarly compare the (to — > 0) dc conductivity ( as well as its normal state oj = inverse, 
a dc (^) = / 9 (^')' tne resistivity) with the counterparts in the high T c superconductors and find reasonable agreement with trends 
in doping and temperature. While condensed matter "simulations" based on these atomic Fermi gases are now of wide interest, 
our approach here is somewhat different because we focus on condensed matter phenomena which currently exist in nature, for 
example the pseudogap in the cuprate superconductors. We have emphasized the commonality for these two systems in past 
reviews 6 and more recently 7 in the context of analogous tools such as radio frequency spectroscopy, which is closely related to 
photoemission 10 . It seems natural then to investigate here the analogies in transport. 

Our work is organized around a second important premise: any theory of transport or scattering in the superfluid Fermi gases 
(both above and below T c ) must be formulated in a fully gauge invariant and sum-rule consistent fashion. One could imagine 
doing more and more sophisticated many body theories than standard BCS-related approaches (i.e. summing more and more 
diagrams) but, if sum rules and conservation laws cannot be demonstrated to hold, we would argue the significance of these 
calculations is problematic. Indeed, a major theme of our recent 1415 work which is summarized in this overview is to arrive at a 
demonstrably consistent theory of this exotic superfluidity by extending the important contributions of Nambu 16 to the case of 
BCS-BEC crossover. This will serve as a basis for computing gauge invariant and properly conserving response functions which 
are measured in transport and scattering experiments. It is not as straightforward to demonstrate this consistency in Monte Carlo 
or other numerically-based schemes. Moreover, as theories of BCS-BEC crossover become more and more diagrammatically 
complicated it is not straightforward to determine that they are compatible with gauge invariance and conservation laws. Once 
this compatibility is established the next step is to assess the physical implications of this consistency in experimental probes. 
We stress here that our approach to BCS-BEC crossover is semi-analytical and therefore reasonably transparent. Moreover, 
unlike alternative (more numerical diagrammatic schemes 17 ) the normal to superfluid transition is properly smooth and second 
order. 

A systematic treatment of transport in the Fermi superfluids requires the introduction 18 of inter-dependent fermionic as well 
as bosonic (or pair) excitations. In past transport literature there has been a focus on either one 19 or the other 4 , but not both. 
Here we use a Kubo-based formalism which readily accommodates the simultaneous bosonic and fermionic excitations of the 
normal state and condensate and thereby addresses n e f / quite accurately while the alternative Boltzmann or kinetic theory-based 
approaches do not naturally incorporate these multiple statistical effects. The Kubo approach includes scattering processes via 
the lifetimes 20 which appear in the various Green's functions, while Boltzmann schemes treat lifetimes via collision integrals. 
However, because the physics of this dissipation is principally associated with the many body processes of boson-fermion inter- 
conversion (via a parameter 7 which appears throughout this paper), it should be satisfactorily addressed only in theories which 
treat the mixed statistics. 

We stress in this paper that dissipative transport in the normal phase must also be addressed if one is to fully understand the 
counterpart below T c case. Except in strict BCS theory, this normal phase is affected by the below T c pairing correlations. As a 
result, consistency checks need to be applied above T c and the behavior must necessarily reflect that the transition is second order. 
Our central theme is based on the fact that the fermionic excitation gap or pairing gap persists above T c as a "pseudogap" (except 
in strict BCS theory); this leads to important transport implications such as "nearly perfect fluidity" (reflecting anomalously small 
viscosity 1 ) or in the analogous charged system "bad metal" 2 behavior (reflecting anomalously small conductivity). 

To add to the case for simultaneously studying both the viscosity and dc conductivity we note that the wealth of data available 21 
and the relative ease of measurement (compared to those in atomic or RHIC experiments) make the cuprates a particularly useful 
analogue for the nearly perfect fluids. In addition to dc measurements, one can probe 21 the conductivity as a function of 
frequency cj, <r ac (uj), over a wide range of to. Because of the existence of a frequency sum rule, theories of the (dc) conductivity 
are thereby highly constrained. Importantly, this sum rule and an in depth understanding of the conductivity serve to constrain 
analogous microscopic theories of viscosity. Conventionally, "bad" metals are systems in which the estimated mean free path 
I is shorter than all length scales; along with anomalously low conductivity this leads to the absence of resistivity saturation. 
The descriptive "perfect" is also associated with a situation in which I is small compared to physical length scales 22 . In strongly 
correlated superfluids, we re-iterate here that small I does not solely reflect short transport lifetimes t but rather a notable 
suppression in the effective "carrier number". The influence of bad metallicity on superconductivity was studied in Ref. 2. Here 
we emphasize the converse: the influence of superconductivity on transport 23 . Moreover, as will be discussed in more detail 
below, pairing fluctuations and the phase fluctuations invoked earlier 2 are clearly distinct. 



3 




u.o i 1,3 

T/T c T [KJ 



Figure 1: Viscosity in helium-3 (left) from Ref.28 and helium-4 from Ref.29 (right). The solid curves on the left represent theoretical 
predictions which were not verified. The viscosity vanishes for low temperature in the case of helium-3 because of a suppression of fermionic 
excitations. The behavior of r\ in helium-4 is governed by the rotons and phonons. 
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Figure 2: Temperature dependence of lj — > conductivity 30 (left) and resistivity 31 (right) in the cuprates. The conductivity exhibits a maximum 
as a function of temperature because of competition between the effective carrier number and quasiparticle lifetime. Increasing the doping 
tends to lower the resistivity. 



We have included a summary of the main results of this paper in Section I C. For ease in reading this paper, the reader who is 
not interested in the technical details can skip Sections II and III and go directly to Section IV. 

A. Experimental Overview 

To understand viscosity in fermionic superfluids, it is useful to begin with helium-3 which has been successfully described 
using BCS-based approaches, albeit for the p-wave case 24 ' 25 . Here experiments 26 indicate that 77 drops off rapidly to zero in the 
superfluid phase. This is shown in Figure 1 to the left. Interestingly, there is a minimum in 77 above T c which is associated 
with strict Fermi liquid behavior 27 . In a Fermi liquid the number of carriers and mass are both T independent, while the inter- 
fermion scattering lifetime varies as T~ 2 . The standard interpretation of the data below T c is that 77 decreases with decreasing 
T as a result of the suppression of fermionic excitations at low T. In a strong magnetic field one spin component of the triplet 
is driven normal and this leads to a very different behavior for the shear viscosity 28 , in which (even below T c ) it reflects the 
normal Fermi liquid behavior above T c . In this A\ phase, the low temperature behavior exhibits an upturn at low T; this is not 
to be associated with coupling to collective modes or phonons, but rather reflects a residual normal component. In BCS-based 
superfluids, we stress 16 that Nambu-Goldstone boson effects do not naturally enter into the transverse transport properties such 
as 77. By contrast, in the helium-4 counterpart shown on the right in Fig. 1, the single particle bosonic excitations couple to the 
collective (Nambu-Goldstone) modes, leading to an upturn 29 in 77 at low T, which has also been predicted (but not seen) for the 
atomic Fermi superfluids 4 . 

For strongly correlated charged systems, the counterpart experiments are summarized in Figure 2, here for the high T c 
cuprates, which have <i-wave pairing. On the left is the low T, ui — » conductivity for a T c m 9QK sample. This shows 
the fairly generic maximum below T c . The figure on the right is the resistivity or inverse conductivity for different stoichiome- 
tries as a function of temperature. Above T c , two crucial points are that p = a^ 1 is nearly linear with T and its magnitude seems 
to decrease from the UND (underdoped) to the OVR (overdoped) samples. It is inferred 32 from similar systematic studies that 
the effective number of carriers at fixed T is substantially depressed, varying as the doping x rather than the expected 1 + x. 
Importantly, for the present purposes a(T) seems to approach zero at the lowest temperatures. The latter point is consistent with 
the vanishing 77 shown in the previous figure for the fermionic superfluid. The pronounced maximum is thought to arise from the 
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Figure 3: Viscosity experiments on unitary Fermi gases from recent experimental results showing the shear viscosity r\ and the ratio of viscosity 
and entropy density r)/s, Ref. 12. Importantly, the viscosity vanishes are low temperature and this puts a strong constraint on the theories of 
the unitary Fermi gas. 



competition between the decrease in n c g and the increase in the fermionic lifetime at low T. Such a competition is not nearly as 
apparent in an s-wave system where n e g is exponentially suppressed. 

In Fig. 3 we summarize recent viscosity experiments from the Duke group 12 on a unitary trapped Fermi gas. One sees here 
that despite previous predictions, the viscosity and its ratio to entropy density are both strongly suppressed at low T. We will 
return to all of these figures throughout the paper. 



B. Theoretical Overview 



The BCS-BEC crossover theory that is adopted here is based on a natural extension of the BCS ground state, 

*o = n k (u k + ^ k 4 t cL H )|o). (l) 

This ground state is often called the "BCS-Leggett" state. The observation that this state is much more general than originally 
presumed forms the basis for BCS-BEC crossover theory 33 35 . To implement this generalization of the BCS ground state, 
all that is required is that one solve for the two variational parameters u k and w k in concert with a self consistent condition 
on the fermionic chemical potential \i. As the attraction is increased, fj, becomes different from the Fermi energy Ep. The 
two variational parameters w k and u k can be converted to two more physically accessible parameters associated with the zero 
temperature gap (or equivalently order parameter) and fi. We stress that while in the cuprates the system is far from the BEC 
limit (in large part because of tf-wave lattice effects 36 ), it is nevertheless quite distinct from the weak-coupling BCS limit due to 
the anomalously short coherence length 9 . 

Since this wave function is so closely related to the BCS state, it is natural to ask whether its behavior away from T = can 
be consolidated into as simple a formalism and physical picture as there is for the ground state. In Section II, we answer this 
question in the affirmative by recasting the equations of conventional BCS theory using a formalism that can then be readily 
generalized to include BCS-BEC crossover. 



1. Simple Physical Picture of BCS-BEC Crossover Scenario 



Before introducing the T ^ formalism, we present a simple picture of the excitation spectra which ultimately enter into 
transport. The top row in Fig. 4 (from left to right) shows the schematic behavior as one passes from the T = BCS-Leggett 
ground state to the above T* Fermi liquid. The red (dotted circles) pairs are associated with net finite momentum, while the 
blue (solid circles) pairs correspond to the phase coherent condensate with zero center of mass momentum and the lone arrows 
represent fermionic excitations. The first panel shows that the ground state consists of fully condensed pairs [as in Eq. (1)], while 
the second panel shows that below T c but above T = there are both condensed and non-condensed pairs along with fermionic 
excitations. These non-condensed pairs persist above T c (third panel) in the form of "preformed" pairs, while the condensed 
pairs are no longer present. Finally at temperatures above T* all bosonic-like excitations are absent; the only excitations are 
fermionic. The second panel with < T < T c is the most interesting from the perspective of the present paper. In the cuprates, 
this is the regime in which the widely discussed 37 "two-gap" physics appears. Here the coexistence of the condensate and of 
non-condensed pairs, leads to two gap contributions 38 , one associated with the pseudogap (called A pg ) and another associated 
with the condensate (called A sc ). 

In contrast to this pairing fluctuation picture of the BCS-BEC crossover theory stands the phase fluctuation picture shown in 
the bottom row of Fig. 4. In this scenario of the pseudogap 39 there exist finite size regions of superconducting order in the normal 
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Figure 4: (Color online) Cartoon comparing the BCS-BEC crossover and phase fluctuation scenarios. Throughout, blue closed circles, lone 
arrows, and dashed red circles represent condensed fermion pairs, unpaired fermions, and finite lifetime pairs, respectively. The crossover 
theory is distinguished by the presence of noncondensed q / pairs for nonzero temperatures less than T* . The defining feature of the phase 
fluctuation picture is the presence of different phase domains above T c , indicated by the regions labelled with different $'s. 



state. While the amplitude of the pairing gap Ak is fixed in these regions, the superconducting phase, labeled by $, fluctuates 
spatially between them. The degree of phase fluctuation increases with temperature, and the Fermi liquid state is reached at 
T > T*. In contrast, the phase is long-range ordered below T c and the superconducting phase is described by the BCS state 
consisting only of condensed pairs and unpaired fermions. In this picture there appears as yet to be no counterpart to the "two 
gap physics" below T c seen in the crossover scenario. In a simplified fashion, it could be said that the preformed pair theory 
is concerned with fluctuations in momentum space (in terms of q ^ pairing) while the phase fluctuation picture focuses on 
fluctuations of the phase in real space. 

In the BCS-BEC crossover scenario, the effects of condensed [i.e., superconducting (sc)] and non-condensed [i.e., pseudo-gap 
(pg)] pairs are described by two distinct contributions to the fermionic self energy E(k, to) 



S(k,w) = £ 



pg-M 



A 2 . 



A L,k 



(2) 



The gap functions A pg k and A SCi k are assumed to follow either a simple s or d-wave form. The condensed pairs have the usual 
BCS self energy contribution, S sc , while the self energy of the non-condensed pairs S pg possesses an additional term, 7, with 
7 _1 reflecting the finite lifetime of the non-condensed pairs. This form of E pg was derived microscopically in Ref. 40 using a 
T-matrix approach (see below). It plays a central role in transport, largely through Ward identities which relate the self energy 
to transport properties. At the microscopic level, it is important to stress that the above expression for S pg is not generic to 
all T-matrix theories, but strongly depends on an underlying BCS-like structure of the ground state associated with the present 
approach. 



2. Lifetime Effects 



In our discussion of dissipative transport, a crucial point is to address the origin of finite lifetime effects. Throughout this paper 
we will argue that the central dissipation process is associated with the inter-conversion from fermions to pairs. We stress that 
this is a many body effect and should not be associated with the two body scattering length. Nevertheless, at unitarity, these inter- 
conversion processes are likely to lead to the shortest lifetimes simply because the number of fermions and bosons is roughly 
equal there. This is in contrast to the BCS (BEC) regime in which there are virtually no bosons (fermions). This physical picture 
is consistent with Eq. (2) which above T c has been rather widely adopted by the high T c community 41 ^ 3 and the cold Fermi gas 
community 10 . In this second context, it is this form of the fermionic self energy (or equivalently spectral function) which is to 
be associated with the downward dispersing quasi-particles revealed in momentum resolved radio frequency experiments 7 44 . In 
the cuprates, above T c , it is the finite lifetime of the non-condensed pairs which leads to the interesting physics associated with 
the "Fermi arcs". These have been interpreted as a blurring of the d-wave nodes. Importantly, below T c one sees their sudden 
collapse to conventional point nodes 43 45 as a result of the onset of the order parameter A sc . 
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Figure 5: Summary figure of d-wave conductivity (top 4 panels) and s-wave viscosity (bottom 4 panels) as a function of temperature. Both 
panels (a) and the resistivity at the top in panel (d) (or inverse conductivity) can be compared with experiments in Figures 3 and 2. Panels (b) 
are contributions to transport from 3 components. Here red = pg, blue = sc, while green reflects the difference to make up the total. Panels (c) 
plot the effective carrier number and (d) the inverse transport coefficients. 
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Figure 6: (a) Comparison of shear viscosity r\ = anh and experiments 46 (red triangles) at unitarity for a trapped gas. In theory plots (black 
dots) we use the calculated thermodynamics for the trap energy E and entropy density s. The inset in (a) plots the estimated transport 
lifetime from radio frequency experiments 13 , as compared with the experimentally deduced lifetime needed for an exact fit to the theory, (b) 
Comparison of r//s. The blue dashed line labels the quantum lower limit of 77/s given by string theory 3 . The last panel on the right converts 
the horizontal axis in (b) to temperature. 



C. Central Results of this Paper 

1. Summary of Central Formulae 



In the weak dissipation limit (where 7 is small), we may write down a combined expression for the shear viscosity and dc 
conductivity. We have 



o- J C ^J dp m 2 El [Up T Vp) [ dE p ' lha {i) 

Here E p is the usual BCS dispersion, E p = ^/£ p + A 2 where A 2 = A 2 9 + A 2 C , and f(E p ) is the Fermi function. We define 
c v = l/lbn 2 ,^ = 1/3-7T 2 and the "coherence factors" are u 2 , u 2 = |(l±-g). The transport lifetime r, )(T scales inversely with 
7. This reduces to the usual BCS results when A 2 g = 0. Eq. (3) contains two strong correlation effects which reduce the size 
of 77 and a: the presence of a gap A(T) in E, which persists into the normal state and the presence of the factor — A 2 9 which 
reflects the reduction in the number of fermions contributing to dc transport as a result of their conversion to bosons. While the 
gap is fermionic in nature, this term is associated with bosonic degrees of freedom. In the stronger dissipation limit, the shear 
viscosity and conductivity can be rewritten in terms of generalized Green's functions (to be defined more precisely in Sec. III). 
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(4) 



where P + = (i(u> n + f!(),p + P~ = (iu! n ,p — and w + = cj + i(J + . Note that while one can interpret F sc as the 
usual Gor'kov Greens function reflecting superconducting order, there must also be a counterpart F pg (discussed in detail later) 
which reflects non-condensed pairs. In previous work on both cold gases 47 and cuprates 48 ' 49 only the first term involving GG 
was included, and, moreover, this was correctly recognized 47 49 as inadequate. 



2. Summary of Transport Results 

We next present summary numerical plots and compare the viscosity and c?-wave conductivity calculations based on Eqs. (4). 
For definiteness we take a simple Fermi liquid j(T) oc T 2 temperature dependence for both transport properties. An essential 
change from s-wave to d-wave pairing is associated with the fact that for the latter the carrier number is no longer exponentially 
activated. As a consequence, the transport behavior tends to be more metallic. 

The four top panels in Fig. 5 plot (a) the T c -normalized, to — > 0, rf-wave conductivity, Re<r(T/T c ). In Fig. 5 (b) we plot a 
decomposition of this conductivity with three different contributions as color coded, so that the condensate (sc) is blue and the 
pg contributions are red, and what makes up the difference is shown in green.. Plotted in Fig. 5(c) is the effective carrier number 
(n/m(T)) e ff which is obtained by multiplying a by j(T) corresponding to the inverse lifetime. Finally in Figure 5(d) we plot 
the resistivity, or inverse conductivity, above T c . The lower panel of figures corresponds to the counterpart plots for the s-wave 
viscosity, in terms of a defined via rj = anh. 

The figures in the (b) column show that the pg contributions are negative in both cases while the condensate gives a positive 
contribution to a and a negative contribution for rj; this can be directly seen from Eqs. (4). The remainder (associated with the 
addition of the GG terms in Eqs. (4)), which yields the respective totals, is shown in green. That the (pg) contribution from the 
non-condensed pairs (red) lowers the conductivity and viscosity is because the presence of non-condensed pairs means fewer 
fermions are available for dc transport. Plotted in in the (c) column one sees that the effective carrier number (n/m(T)) e g 
increases more or less monotonically as temperature is raised. This represents a generic figure, since here the temperature 
dependence of the lifetime has essentially been removed. Thus, this figure shows that the effective carrier number associated 
with the conductivity and its shear viscosity counterpart are increasing functions of temperature, strongly suppressed by both the 
presence of afermionic gap and the presence ofbosonic degrees of freedom. In Fig.5(d) we plot the inverse transport coefficients. 

The first two figures (top and bottom) on the left in Fig. 5 seem to capture the qualitative experimental transport features shown 
in Figures 2 for the cuprates and Figure 3 for the cold gases. The conductivity exhibits 30 a maximum below T c , as can be seen 
experimentally from Fig. 2, while the viscosity coefficient monotonically increases, as observed 12 . As shown in panel (d), the 
conductivity of the normal state is appropriately metallic, but suppressed by the excitation gap. Importantly, the resistivity has a 
nearly linear temperature dependence as observed experimentally, in Fig. 2. 

For the cold gases, more direct comparison with experiment involves inclusion of the trap. As is conventional 6 , we include 
trap effects via the local density approximation (LDA), or, equivalently Thomas-Fermi approximation. Figure 6 presents a 
comparison of the viscosity coefficient a between theory (based on the RF-deduced lifetime), as black dots, and experiment 46 
(red triangles) as a function of E. These calculations can also be compared with more recent experiments summarized in Fig. 3. 
Our calculations, which predated the latest experimental data, will be discussed in more detail later. One can see, however, that 
the trends are compatible. 



II. THEORETICAL FRAMEWORK OF BCS-BEC CROSSOVER 



A. T-matrix Derivation of BCS Theory 



In order to understand how to address BCS-BEC crossover at finite T, we now rederive standard BCS theory from a T-matrix 
scheme. Important here is that BCS theory can be viewed as incorporating virtual non-condensed pairs. Here we consider the 
general case applicable to both s and d-wave pairing by incorporating ip^ = [cosf^) — cos(k y )] for the latter and taking it to be 
unity for the former. These virtual Q ^ pairs are associated with an effective propagator or t-matrix which is constrained to 
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be of the form 



<(Q) = i+rrv r r ^ • (5) 

1 + U l^K G ^ G 0-Kr+Q</? k _ q/2 

in order to yield the standard BCS equations. This t-matrix is associated with a summation of ladder diagrams in the particle- 
particle channel [see Fig. 7 (a)] and importantly depends on both G and Go, which represent dressed and non-interacting Green's 
functions, respectively. We use K, Q to denote four- vectors with K = (iui n , k), Q = (ifli, q), and where ui n and 17; are fermion 
and boson Matsubara frequencies, respectively. In order to describe pairing in the d x 2_ y 2-wave channel, we write the attractive 
fermion-fermion interaction in the form E/k.k/ = t^k^k', where U is the strength of the pairing interaction. As in bosonic 
theories, non-condensed pair excitations of the condensate are necessarily gapless below T c . This means that t(Q — > 0) — > oc 
which is equivalent to the vanishing of the effective pair chemical potential, fj, pa i r = 0, for T < T c . Thus we have a central 
constraint on the T-matrix 

t-\Q -> 0) = -> n pair = 0, T < T c (6) 

In order to show that the above condition is identical to the BCS gap equation, we need the appropriate form for Gk- In the BCS 
theory the fermionic self energy that appears in the fully dressed Green's function, Gk, is of the form 

^sc,K = ^t sc (Q)Go-K+Q¥l-c i /2 = -^A 2 sck S(Q)Go-K+Q = -A 2 s C MGo-K (7) 

Q Q 

where A SCj k(T) = A sc (T)ipk is the superconducting order parameter. The full Green's function is then obtained via the Dyson 
equation, Gk = [Gqk ~ ^sc,k] _1 , which, when inserted in Eq. (5) yields the BCS gap equation below T c 

1 = -"E lz iPrf < 8 > 



with E^ c = + A? sc k . We have, thus, used Eqs. (5) and Eq. (6) to derive the standard BCS gap equation within a T-matrix 

language and the result appears in Eq. (8). Equation (6) above can be viewed as representing an extended version of the Thouless 
criterion of strict BCS which applies for all T < T c . This derivation leads us to reaffirm the well known result 50-52 that BCS 
theory is associated with one bare and one dressed Green's function in the pair propagator. 



B. Generalization to BCS-BEC crossover 



To address BCS-BEC crossover, we presume that the non-condensed (Q ^ 0) pairs are no longer virtual. Thus the T-matrix 
of Eq. (5) in general possesses two contributions: the q = contribution that gives rise to the formation of the condensed or 
superconducting pairs and the q ^ contribution that describes the correlations associated with the non-condensed pairs. As a 
result, the fermionic self-energy also possesses two contributions which are given by 



^( K ) = X! t (Q) G 0,-K+Q l fi'l- cl /2 = y^sc(Q) + t pg(Q)}Go,-K+Q l pl-q/2 = ^sc.K + ^pg,K (9) 

Q Q 

The resulting full Green's function, G _1 = Gq 1 — T, sc — T, pg is illustrated in Fig. 7(b). While as before, E sc ,k = — k Go,-if, 
we find numerically 53 ' 54 that ^ P9 .k is in general of the form 

A 2 . 

Z pg . K « - j** . (10) 
to + Ck + 11 

with A ra ,k = Apgtpk. That is, the self-energy associated with the non-condensed pairs possesses the same structure as its BCS 
counterparts, albeit with a finite lifetime, 7 _1 . Physically this arises from the fact that t pg (Q) is strongly peaked around Q = 
below T c where the pair chemical potential is zero and for a range of temperatures above T c as well where this chemical potential 
is small. 

Analytic self consistent equations for A pg and A sc can be obtained microscopically when we consider the small 7 limit where 

£(#) « -(A 2 c , k + A 2 9 , k )G ,-K = -A 2 G ,-k (11) 
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with 

^ s = -Ewo) (12) 

Q 

Here, we have used the fact that because of the vanishing of [i pa i T below T c , the bulk of the contribution to E ps in the ordered 
state comes from small Q. This then leads to an effective pairing gap A(T) whose square is associated with the sum of the 
squares of the condensed and non-condensed contributions 

Al(T)=Al ik (T) + Al gM (T) 

Note that the full gap Ak remains relatively T-independent, even below T c , as observed, because of the conversion of non- 
condensed (Ap g .k) to condensed (A SCj k) pairs as the temperature is lowered. 

The gap equation for this pairing gap, Ak(T) = A(T)<^k, is again obtained from the condition t~ g {Q = 0) = 0, and given 

by 

This analysis can be made more explicit after analytical continuation so that 40 , t pg (u), q) ~ [Z(il — f2 q + \x pa ir) + ^qI ^ 
where Z = (d x /dn)\n= , q =o, ^ q ~ q 2 /(2M h ) with the effective pair mass M^ 1 = (l/2Z)(d 2 x/dq 2 )\n= , q =o- Near T c , 
Tq — > faster than q 2 as q — > and will be neglected. Then A 2 g « Z^ 1 6(O q — /J, pa i r ). 

Note that one needs to self-consistently determine the chemical potential, /j,, by conserving the number of particles, n = 
2 J2k Gk, which leads to 



K 



1 - |£- + 2^/(£ k ) 

^k 



(14) 



Eqs. (12), (13), and (14) present a closed set of equations for the chemical potential [i, the pairing gap Ak(T) = A(T)<^k, 
the pseudogap A pgj k(T) = A pg (T)ip k , and the superconducting order parameter A sc ,k(T) = A sc ipk with A SC (T) = 



A 2 (T) - A pg (T). Following this approximation, A pg {T) essentially vanishes in the ground state where A = A sc . This 

is to be expected from the BCS-Leggett wavefunction in Eq. (1). In this way, the "two gap" physics disappears in the ground 
state. Importantly, numerical studies 36 show that for rf-wave pairing, there is no superfluid phase in the bosonic regime where /i 
is negative; the pseudogap is, thus, associated with the fermionic regime. 

III. TRANSPORT THEORY AND GAUGE INVARIANT APPROACHES TO SUPERCONDUCTIVITY AND SUPERFLUIDITY 

Our transport theory for BCS-BEC crossover is based on linear response theory for both the density (or charge, labelled 
C) and spin (labelled S) degrees of freedom. In this approach, the U(l) electromagnetic (EM) gauge symmetry and the spin 
rotational symmetry around the z axis play important roles and an understanding of transport in strongly correlated superfluids 
has to incorporate in a central way the related conservation constraints. These enter via (i) the transverse f-sum rule. Application 
of the latter to the conductivity is, in turn, related to (ii) the absence (above T c ) and presence (below T c ) of a Meissner effect. 

The perturbing Hamiltonian can be written in a compact form J d 3 r(\'^'^,tp'l'ip a > + h.c.) where \ a ' a , oc 5 aa > and A^ CT , oc 
5 a a'- Here ip^itp) are the fermionic creation (annihilation) operators, a =t or |, f= — 4- an d & = We represent the 
density- density, current-current and spin correlation functions as, Xpp< X.JJ an d Xss- Experimentally, the last of these three 
can be probed by spin-preserving and spin flip (two photon) Bragg scattering. The shear (77) and bulk viscosities ((2) and the 
conductivity may be written in terms of the longitudinal (xl) an d transverse (xt) components of the current-current correlation 
functions x.u 

77 = — m 2 lim lim -— Imyrfw, q), (15) 



4 
—7 

3 ' w-K) q^o q 2 



(2 + 7;V = -m 2 lim lim — Imxi(q, w), (16) 



*{ U 0) = - lim lim ImXT(u; ' q) (17) 

w->0 9 ->0 0J 
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where the longitudinal \l = q • XJJ 1 <3 and transverse \t — {J2a=x Xjj ~ Xl)/2 susceptibilities satisfy 8 



Um p^ Im^q) = njT)^ 

+oc dm ^ ImxL(^,q) ^ = n_ ^ 



to m 



f 

J — ( 



dflReaiQ) = e 2 — (20) 
m 



Because we simultaneously discuss both neutral and charged systems, it is useful to define e = 1 for the neutral case. Here n n 
is the particle number of the normal component in the superfiuid and the sum rules are explicitly written for the Fermi gas in 
which there are no bandstructure effects. For scattering probes, we define the associated structure factors for spin and charge or 
density in terms of closely related response functions, x pp and xss'- 

5 c (q,w) = coth(^)Imx pp (q,cj) and 5 s (q,w) = coth(^)Imxss(q, w) 

TT 11 IT 11 

Since the conservation laws for particle number and spin, d^J pi = and <9 M J^ 5 = are satisfied, the following two sum rules 
must be respected at all temperatures in the whole BCS-BEC crossover regime: 



f 

j — i 



n c^ 2 

dojujSc s(q, w) = with lim 5cs(q, w) = 0. (21) 



The EM kernel is defined by J = —K ■ A, where K{Q) = e 2 (n/rn) + P{Q), and the paramagnetic contribution, given 

by P(Q), is associated with the normal current resulting from fermionic and bosonic excitations 14 ' 55 . In the superfiuid phase, 
the density correlations which enter into scattering and the current correlations which enter into transport can be schematically 
written as a sum of 3 terms where for convenience we drop the (u>, q) arguments: x.JJ = +tolij and \pp = P 00 +Coll p 

The counterparts for the spin degrees of freedom are xss — Qs + — an d Xss — Qs°- Here P, P 00 and Qs, Qg represent the 
"bare" contributions. Collective mode effects in the charge response, which are not present in the spin response, must also be 
included in the longitudinal response below T c . These appear in the above equations as Soli j and Collp. These collective mode 
effects are essential for insuring that the sum rules and related conservation laws are satisfied. 

In the most general case, the diamagnetic current is expressed in terms of the inverse band mass <9 2 £k/ 'dk a dkp (with a, /3 = 
x,y,z), via (n/m)di a = 2 ^2k a ( < 5 2 £k/<9k<?k)G.R". Importantly, the latter contribution which is temperature independent, 
should not be confused with (n/m(T)) e g. We stress that this effective carrier number is sensitive to the pairing gap A(T), while 
the diamagnetic contribution is not. We integrate the expression for the diamagnetic contribution by parts and use the self energy 
equation and the Generalized Ward identity to obtain (See Appendix A) an alternate form above T c 



, m / dia ^ <9k <9k 



(22) 



This exact t-matrix based equation is significant because it has cast the diamagnetic response in the form of a two particle 
response function. That there is no Meissner effect in the normal state is related to a precise cancellation between the diamagnetic 
and paramagnetic terms. Noting P(0) = — e 2 (n / m) di a , we may extend P(0) to finite Q to infer 



P(Q)=2e 2 J2- 

r)U r)k \GkGk+Q + *^2t pg (P)Go t P-K-QGo t P-KGK+QGK 



dk dk 



(23) 



Thus far, our discussion has been quite general, and we have circumvented any discussion of specific transport diagrams by 
building in the absence of a Meissner effect above T c . One can alternatively 14,55,56 introduce the Aslamazov-Larkin (AL) and 
Maki-Thompson (MT) diagrams to arrive at the above equation (see Fig 7), but the former which involves two factors of t pg , at 
first sight, appears more complicated. 

Collective mode effects are not present in the viscosity and conductivity, because they represent transverse probes. However, 
they play an important role in the density-density response. In this regard, it is convenient to define 



/•OO 

S± = (Sc ± S s )/2 from which it follows that / dujujS- (w, q) = 0. 

Jo 



(24) 



11 




Figure 7: Characteristic diagrams which enter into the spin and charge response kernels (left). Heavy and light lines correspond to dressed and 
bare Green's functions. Here the wavy lines represent the pair propagator t. The corrections to the vertex include the Maki Thompson (MT) 
and Aslamazov-Larkin (AL) diagrams. Temperature dependence of S_(q) (see Eq. (24) in units of ti/2Ef) for fixed momentum transfer 
q = 1.5&F and frequency u — 2.0Ef at unitarity (right) . This figure shows that S- behaves like an order parameter 14 . 



This latter is a very unusual sum rule. However, it must be satisfied in any consistent theory of superfluidity, providing spin 
and charge are conserved. In Appendix B we show how this sum rule is satisfied above T c for an alternate BCS-BEC crossover 
theory introduced by Nozieres and Schmitt-Rink (NSR) 35 . 

There are other physical consequences which can be deduced once one has a conservation-law consistent theory. The conduc- 
tivity and the shear viscosity can alternatively be written in terms of the bare response P{Q) so that 

cr(w) = -lim with v = -to 2 lim lim — lmP xx (w, q). (25) 

Moreover, in this way, one can see how closely related they are. We will not in this paper address the bulk viscosity, principally 
because we do not have the same level of theoretical control to satisfy the longitudinal /-sum rule, which is also more problematic 
below T c , where the Goldstone bosons appear 16 . 



A. Weak Dissipation Regime 



In the previous section we have set up the general linearized response theory. We now discuss the detailed implementation of 
this formalism based on Eq. (2). Consistency with conservation laws requires that transport diagrams assume a specific form, 
reflecting the behavior of the self energy. We begin with the so-called "weak dissipation" regime, where in Eq. (2) we presume 
the quantity 7 « 0. The strong dissipation regime is discussed in the next subsection. We have found that for the viscosity there 
is very little difference between weak and strong dissipation-based calculations, but this is not true for the d-wave conductivity, 
because of the easy excitations of fermions around the nodes. We indicate in Figure 7 the characteristic diagrams which enter 
into the generalized correlation functions. As a result one can show 15 that the "bare" contributions to generalized charge and 
spin susceptibilities take on an extremely simple form 



P 0C V,q) 

g° s °( W) q) 



E 



\ E $ 



E~ E+E~ 
p p p 



£+£" ± A 2 

Sp Sp - 1 - "si 



A 



pg 



L E+E p 

E i 



uj 2 -(E++Ep) 2 



(l-/(25+)-/(25-)) 



Ep E+E p 



£+£- =c A 2 + A' 

>p >p T "sc 1 ] 



IV 



E+E p 



to 



2 -(£+-£ p -) 2 



(f(E+)-f(E p )) 



(26) 



£(w,q) 
<3s(w,q) 



\ - PP 

" TO TO 



£p + - 



=F A 2 

Sp Sp \ ^sc 



A 2 
pg 1 



E + Ep 
E t 



l-f{E+)-f(E~)) 



E p E t E P 



E+Ep 



(27) 



where ui implicitly has a small imaginary part and we define the quantities E+ = E v ± q / 2 and ^ = £ p ± q /2- Note that Eqs. (26) 
and (27) are exactly the same as their BCS counterparts when A 2 S — > 0. It is quite remarkable that when we compare these two 
equations we see that the charge response functions reflect the difference A 2 C — A 2 fl while the spin response functions depend 
on the total pairing gap A 2 = A 2 C + A 2 ff . We may say that spin correlation functions know only about pairing whereas their 
charge counterparts reflect also coherent superfluid order. 
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^From the definition of Xt(w, q) the shear viscosity is given by 



P 2 '<« 2 Sl,, „„^E+E--(+(- -JA 2 



p 



p+ p- I t+ 1— I X A 2 

-(/(£+) - /(£?-)) p % + X P (W - 



2£+£ p 



(28) 



where we have used the abbreviated notation: 5i(o;) = — — E~ ), ^(w) = <5(w — 1?+ + S~ ). Physically, two types 
of terms appear in the above equation. Both are well known in standard BCS theory. The first refers to processes which require 
a minimal frequency of the order of 2A(T); they arise from the contribution of fermions which are effectively liberated by the 
breaking of pairs. The second of these terms, involving 82, arises from scattering of fermionic quasi-particles and is the only 
surviving contribution to the viscosities, which are defined in the lo — > limit. Note that both contributions involve the difference 
of the condensed and non-condensed components (SA 2 = A 2 C — A 2 g ) with opposite overall signs. The low to quasi-particle 
scattering processes are reduced by the presence of non-condensed pairs - because they are associated with a reduction in the 
number of fermions. By contrast in the high w w 2A limit the number of contributing fermions will be increased by breaking 
pairs. We next take the low u>, q limits in Eq. (28) and introduce lifetime (or dissipation) effects by writing delta functions as 
Lorentzians S(uj ± q • V_E P ) — > - (^ ±q . v 7 ff y^+^ ■ ^ n tn i s wav we arrive at Eq. (3) which was presented in the summary section. 
Note that 77 assumes a form similar to a stress tensor- stress tensor correlation function. 

The transport expressions summarized earlier in Eq. (3) correspond to the generally familiar BCS form 24 25 except for the 

A 2 

effects associated with non-zero A pg which appears as a prefactor 1 — -gff. This deviation from unity can be traced to the AL 
diagrams and hence may be interpreted as extra "bosonic" contributions which contribute to the normal fluid n n . Their presence 
is physically required for the sum rule consistency of \t- Note that these terms enter with an overall reduction coherence-effect- 
prefactor associated with £ 2 /E 2 . The negative sign for this bosonic term comes from the fact that only fermions contribute 
directly to transport within a BCS-Leggett-based scheme; the more pairs which are present, the lower the fermionic contribution 
to the viscosity. We see from Eq. (4) that a and 77 are alike except for momentum power law prefactors and coherence factors, 
When extended to finite frequency we note that unlike in BCS theory, there is a pair creation and pair-breaking contribution to 
the optical conductivity. This term is necessarily absent in a clean BCS superconductor, but may well be the origin of the widely 
observed 21 "mid-infrared" contribution to the finite co conductivity. 

B. Explicit Proof of the transverse f Sum Rule 

An important check on our microscopic scheme is to show that it satisfies the sum rule for the transverse component. The sum 
rule which are going to prove is in Eq. (18). Our proof here is explicitly for the weak dissipation form of the response functions 
and for the Fermi gas, where there are no bandstructure effects. We generalize this later in the paper. 

The total particle number is 

„ . ?(l .i (1 . vw) ).2 ? ^(lz^ A .. 4 ^)). 

The superfluid density at general temperature 38 is given by 

1 2 J2 



n s 



Therefore 



n„ =n-n s 



2 V P 2 ( A gg 1 ~ 2f(E p ) El A 2 g df{ E v ) \ 
El \ to 2£ p to dE p J' 



Using Xjj = P + leads to 

f + °° dco Im XT 2 v p 2 / A p 2 g 1 - 2f(E p ) E 2 p A p 2 g d f(E p ) s n n 
q^aj_ oo it { lu * 3m^El\m 2E p to dE p J to' 
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By contrast, (except in special cases, such as q — > 0, above T c ), the longitudinal sum rule requires numerical proof. Following 
the analysis in Ref. 14 we find agreement with the sum rule to within 5-10%. In addition to sum rule consistency, the appropriate 
diagram set for computing all transport properties must be chosen so that n s (T) vanishes at and above T c . This is somewhat 
more complicated 57 to ensure than in strict BCS theory because there is a finite excitation gap at the transition. This gap or 
pseudogap, in turn, reflects the fact that there are bosonic excitations in addition to the fermions which deplete the condensate. 
In Section III E we present a more general argument for establishing the conductivity sum rule whether the dissipation is weak 
or strong. 



C. Two Photon Bragg Experiments 



In a recent paper we have used this theory to address the dynamical structure factor and thereby show how two photon Bragg 
scattering can be used to establish in situ the presence of coherent order in a superfluid, at any temperature, wavevector and 
frequency. For the most part experiments on unitary gases have relied on sweeps to the BEC, to find evidence for condensation. 
Our analysis is based on the definitions and sum rule in Eq. (24). It, thus, depends on imposing the current conservation laws, 
which have been extensively studied and verified 14 . Using the characteristic diagrams shown in Figure 7, which enter into the 
density-density correlation functions we are now led to an important observation: the quantity S-(oj, q) for all q, to, 1/kpa 
can be used as an indication of in- situ superfluid order, without requiring sweeps to the BEC. We show a plot of this behavior 
in the right-hand panel of Figure 7, where it can be seen that the difference structure factor vanishes in the normal state. An 
interpretation of this figure is that despite previous claims 58 there is no spin-charge separation in the normal or pseudogap state 
associated with BCS-BEC crossover. Spin-charge separation is, however, to be found in the superfluid phase. The counterpart 
normal state calculations can be shown to be valid in an alternate BCS-BEC crossover theory, based on the Nozieres Schmitt- 
Rink scheme 35 ; this is presented in Appendix B. 



D. Strong Dissipation Approach 



We now use the full expression for the self energy Eq. (2) to obtain compatible expressions for transport coefficients The full 
Green's function is given by 

Gk = (iu n - & + i7 - ■ A ^\. -^iV) _1 - (33) 



iUn + £k + H ibJ n + £ k > 

where we have added an extra constant term i-f in order to be consistent with the weak dissipation limit in the case that 7 becomes 
small. To extend P(Q) which appears in Eq. (23) below T c , within a BCS-like formulation one needs to include terms of the 
form F sc kF sc k+q which represent the usual Gor'kov functions as a product of one dressed and one bare Green's function 

(GG ) 

rp A SC .k 1 



iu n + £k 



(34) 



lUJn 



Here, as before, A k = A pg k + A^ c k . Then, in the same spirit as our derivation of Eq. (10) we exploit the fact that t pg {P) is 
strongly peaked at small P, which leads us to approximate 

<9£k+q/2 <9£ k+q / 2 



P{Q) 



K 



where F p „,k 



9k 9k 



GkGk+Q + F SC ^ K F SC ^K+Q — F pg ^KF pg ^K+Q 

G K 



' P9 '" ioJ„ + £k + il 

From Eq. (35) and Re a para (Q) = -Im P XX (Q)/Q the paramagnetic contribution to the dc conductivity 
ReaP ara (0) ^AhnlmJ2[^ (^) 2 (GkG k+q- F pg , K F pg ^ + F sc , K F sc ^) 



^0+ 



In order to be consistent we rewrite Eq. (22), also adding in the usual BCS condensate terms 

W " 2 E(||) [ GkGk +Q F sc ,kF sc , k+ q F pgiK F P9!K+Q 



(35) 
(36) 

(37) 
(38) 



We can use Eq. (35) to arrive at Eqs. (4) which were presented earlier in the form of a summary. Note that in previous work in 
the literature 47 "* 9 only the first term involving GG was included, which was recognized 47 ' 49 as inadequate. 
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Figure 8: (Color online) Theory plots of shear viscosities as functions of T from BCS to BEC, divided by the counterpart unpaired state 
viscosities r] n . The red dashed lines are results in absence of explicit bosonic degrees of freedom (associated with the AL diagrams). In (a), 
T c = 0.127>, T* = 0.13T F and l/k F a = -1. In (b), T c = 0.26T F , T* = 0.507> and l/k F a = 0. In (c), T c = 0.2L7V, T* = 1.287V 
and l/fci?a = 1. Expanded plots in different T regimes are shown in the various insets. 



E. Proof of Conductivity Sum Rule 

We now revisit the issue of compatibility with the important conductivity or transverse sum rule 

/oo 
drtRsa(n) = e 2 (n xx /m) d . a (39) 
-oo 

in a more general fashion. Here we now include bandstructure effects through the effective mass. Note that we must have two 
contributions to the conductivity corresponding to the paramagnetic and diamagnetic terms 

Rea(n) = -ImP xx (Q)/n + n5(n)[ReP xx (Q) + e 2 (n xx /m) d J (40) 

Integrating the first term over frequency we find = — irReP xx (0), while the second (delta function) term yields a term 
+TrReP xx (0), which leaves only the diamagnetic contribution and yields the desired sum rule. Note that this analysis holds 
both below and above T c and that this sum rule is intimately connected to the absence ( above T c ) and the presence (below T c ) of 
a Meissner effect. Importantly, since ( ) dia can be viewed as essentially independent of temperature, when there are approxi- 
mations in evaluating the transport diagrams, it is appropriate to evaluate the chemical potential fi based on the T-independence 
in Eq. (38). It should be re-iterated this this diamagnetic contribution is to be distinguished from n e f / (T), which enters into the 
dc transport. 



IV. RESULTS FOR VISCOSITY 



We now summarize some of our numerical calculations, beginning with the shear viscosity. In Figure. 8 we plot the viscosity 
divided by a "normal state" value r\ n as a function of temperature for 1/kpa = ±1 and for unitarity. This normal state viscosity 

is to be associated with the (temperature dependent) viscosity of the unpaired state = J~ d Pj ^ ( - % \Jr v . In this 

way we take the same lifetimes in numerator and denominator for the plots, which necessarily cancel out of the ratio. For all 
kpa, the ratio drops to zero at low temperatures reflecting the decrease in the number of condensate excitations. The red dashed 
line indicates the behavior when the bosonic excitations are removed and the black solid line shows the full calculated viscosity. 
Thus, the shaded regions correspond to the contribution from the AL diagrams. This contribution is seen to be largest at unitarity, 
as expanded in the inset to Figure. 8(b). The ratio r](T c )/r] n varies from 0.7 to 4.4 x 10~ 5 as one passes from BCS to BEC in 
these three cases. At unitarity the viscosity at T c is reduced by a factor of 10, relative to the unpaired fluid. A weak signature of 
the transition is largest in the BCS regime, as shown in the inset to Figure. 8(a). 

To incorporate trap effects, we begin by summarizing past work on the thermodynamic al properties of cold Fermi gases. 
Counterpart thermodynamic al experiments have played a role in characterizing the viscosity 1 ' 46 ' 59 . Figure 9 presents a compari- 
son of our previous thermodynamical calculations 5 60 with experiments for both the trapped (left) and homogeneous (right) case. 
All data is shown in black, while our theoretical results are shown in red. Agreement is reasonably satisfactory in both cases, 
particularly with the addition of a very small Hartree adjustment which replaces the dashed lines with the solid curves. The 
first four panels (top and bottom, from the left) correspond to comparisons with the trapped case 59 . The next three are for the 
homogeneous case with the black lines corresponding to experiments from Japan 61 . The third from the left two panels compare 
the present theory (in red) with these experiments 61 in the upper panel and other analytical (light blue) 62 as well as Monte Carlo 63 
(dark blue) calculations in the lower panel. The large scale figure shows that the present approach (red curves) is in quite good 
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Figure 9: Left: Calculated 5,60,64 thermodynamical behavior (plotted as red curves, including Hartree shift (solid) and without (dashed). Results 
for the trapped case are shown by 4 left panels and for homogeneous case by 3 right panels. Experiments (black lines) on left correspond to 
those in Ref. 59 and those on right with Ref. 61. Third (lower) panel from left compares with other theories, one of which shows the usual 
spurious first order transition 65 and the other represents Monte Carlo simulations 63 . Here Efg is defined as in Ref. 65 
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Figure 10: Trap profile for viscosity. The various curves show the shear viscosity as a function of trap position for four different indicated 
temperatures and Rtf is the Thomas-Fermi radius. 



agreement with the homogeneous experiments 61 (in black) over a rather wide temperature range. Absent here is the first order 
transition seen in all other analytic theories. 

In Fig. 10 we show the the LDA-calculated viscosity as a function of position in the trap and for various temperatures T. The 
arrows indicate the edge of the superfluid (SF) core. It can be seen that the viscosity is suppressed in the core region. In this way 
the nearly-free fermions at the trap edge dominate the trap averaged value for -q. As a result the calculated viscosity in a trap is 
considerably higher than in the homogeneous case. 

With these thermodynamical calculations and trap profiles as a backdrop, we now return to an earlier figure, Figure. 6, which 
addressed comparison of theory and experiment for the shear viscosity. For a unitary Fermi gas, we estimate the transport 
lifetime from the characteristic broadening of the single particle fermionic spectral function. The latter, in turn, can be found 
here by fitting 7 ' 17 Radio Frequency (RF) cold gas experiments. Recall that this fermionic lifetime is associated with a damping 
contribution in an otherwise BCS-like self energy (conventionally parameterized 40,41 by "7 "). The inset to Figure. 6(a) presents 
a plot of this RF-deduced lifetime as black circles. The red triangles in this inset represent a plot of the lifetime which one would 
infer from the data for precise agreement between theory and experiment. 

Figure. 6(a) presents a comparison of the viscosity coefficient a between theory (based on the RF-deduced lifetime), as black 
dots, and experiment 46 (red triangles) as a function of E. Figure. 6 (b) shows the comparison of r/fs where s is the entropy 
density. We find that r//s appears to be relatively T independent at the lower temperatures. The last figure on the right represents 
a transcription of the horizontal axis in Figure. 6(b) which plots rj/s as a function of temperature, rather than energy E. This 
is based on using the calculated trap thermodynamics 5 to rescale the various axes. Moreover, our calculations incorporate the 
same trap averaging procedure as in Ref. 46. One can anticipate that, particularly at the lower T, the trap-integrated viscosity 
will be artificially higher than for the homogeneous case, since 77 will be dominated by unpaired fermions at the trap edge. It 
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Figure 11: Theoretical plots of resistivity as a function of temperature for three different values of x, corresponding to the pairing onset 
temperature (or strength of attractive interaction). The UND, OPT and OVR labels as the same as those used in Figure 2, from Reference 
3 1 . The resistivity decreases as doping increases, deriving from an increase in the effective carrier number. That the plots stop at temperature 
T = T* reflect the fact that Eq. (2) is no longer appropriate when the system is this far from condensation. 



should be noted that in the homogeneous case the ratio of the viscosity to its normal counterpart is exponentially activated with 
T. This is to be contrasted with the behavior of the entropy which reflects bosonic power laws 5 in T. Overall, it can be seen that 
our calculations agree favorably with the experimental data. Interestingly, the observed behavior appears more consistent with 
previous helium-3 experiments 26 than those in helium-4 29 , as can be seen from Figure 1 . 



V. CONDUCTIVITY IN THE CUPRATES 



It is the goal of this section to address the dc conductivity a(T) for a pseudogapped superconductor such as the high T c 
cuprates. We approach this problem in a fashion which is analogous to our above discussion of viscosity for an ultracold Fermi 
gas. In the cuprate literature one associates a fixed stoichiometry (hole doping) with a given excitation gap at T c , say. The lower 
the hole doping (x) the larger this gap. Since the parameter x is of no particular interest here, we parameterize instead a given 
stoichiometry by the size of A(T C ), or alternatively the size of the temperature T* at which A(T) first vanishes. The two key 
puzzles of the dc conductivity in the cuprates are the near linearity of the resistivity with temperature and the fact that only the 
doped holes (a(T c ) oc x) appear to contribute to transport. 

A central conclusion of our conductivity study is that, just as for the viscosity, the reduction in the effective carrier number 
(n/m(T)) e ff is revealed to play an important role, both in the T and x dependence of transport. This reduction, in turn, is a 
consequence of the presence of an excitation gap which persists into the normal phase, and which increases as x decreases. 
The suppression in the carrier number is substantial relative to the full diamagnetic or sum rule value in Eq. (39). Moreover, 
(n/m(T)) e ff rises nearly monotonically with temperature until T*. This contribution leads to a non-metallic tendency with a 
increasing with T, above T c . In order to yield a metallic resistivity (which increases with T) the contribution of (n/m(T)) e g 
must be off-set by taking ~f(T) to be a higher power than linear. Here we illustrate our results for T > T c with a Fermi liquid 
like behavior ^(T/T c ) 2 , which could plausibly be associated with Fermi arcs 45 , which are extended gapless regions due to the 
smearing out of the <i-wave nodes above T c . This is the most conventional T dependence for transport processes which involve 
inter-fermion scattering. 

We return to Fig. 5 (top panel) which was presented earlier and which summarizes the general behavior. Of particular interest 
is the behavior of the calculated resistivity which is illustrated in Fig.5d which, itself is roughly linear, because of the assumed 
quadratic or Fermi liquid T 2 dependence in 7. Note that these conductivity calculations are specific to the <i-wave case and for 
an s-wave counterpart, it would be very difficult to find metallic behavior. This can be seen by noting that the inverse viscosity 
shown in this figure decreases with increasing T, reflecting the even more strongly suppressed carrier number. 
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Figure 1 1 presents a more detailed plot of these normal state resistivities for dopings that interpolate between the heavily 
underdoped (large T*) and overdoped cases (T* w T c ). The resistivities are normalized by the value p D — p(T c ) for the case of 
highest doping. We may characterize each curve, in order of increasing doping, by the ratio j/A(T c ) — 0.07, 0.20, and 0.35. 
That the size of the resistivities decreases as one increases the hole concentration, largely reflects the change in gap size and thus 
in effective carrier number. Indeed, to a first approximation (n/m(T c )) e g scales close to linearly with x. 

Although it is subtle, one can see from the figure that there is a change in the nearly linear slope with increased doping from 
concave to convex bowing, which may be seen experimentally 31 in Figure 2. In this way, it would appear that the so-called 
"optimal" doping (where the bowing is minimal) represents a form of mid-way point, rather than a specific form of "strange 
metal" 66 . 



VI. COMPARISON WITH LITERATURE 
A. Pseudogap theories 

Since the present transport calculations emphasize the role of the pseudogap, it is useful to clarify the nature of the pseudogap 
which is specifically associated with BCS-BEC crossover. One of the earliest observations that a pseudogap may be present in 
BCS-BEC crossover theories is due to Randeria and co-workers 67 . This first generation analysis focused on "spin gap" effects, 
so that there were a number of claims which are inconsistent with current understanding. Despite statements to the contrary 67 , 
(i) The pseudogap is not a "spin gap", but is also to be associated, as we have seen in the present transport calculations, with 
the "charge channel". This mirrors the experimental observations 68 on the nature of the cuprate pseudogap: "It is a quasiparticle 
gap and not just a spin gap." (ii) The pseudogap is not associated with "spin-charge" separation above T c . In this paper (Section 
III E and Appendix B) we have shown that the absence of a Meissner effect leads to the non-separation of spin-charge degrees 
of freedom, (iii) The pseudogap in crossover theories is quite distinct from the pseudogap associated with the phase fluctuation 
scenario (as discussed in Section IB 1). 

Our group was one of the earliest (i) To flesh out an understanding of the nature of this pseudogap, and, in particular to show 69 
that the BCS-BEC crossover-pseudogap is not only associated with spin as was claimed 67 , but that more generally it represents a 
quasi-particle gap. This observation was based on calculations of the fermionic spectral function, (ii) To address gauge invariant 
electrodynamics in the presence of a pseudogap 38 56 . (iii) To introduce the concept of a pseudogap into the cold gas literature 70 
and (iv) To introduce 40 the now widely used expression 41 for the self energy (Eq. (2)). 

B. Cold Fermi Gases 

There has been a significant effort devoted to reaching a theoretical understanding of the various viscosities associated with 
highly correlated fermionic superfluids 3 . A detailed overview exploring the relation between cold atomic gases and hot quark 
gluon plasmas was presented by Schaefer and Teaney 71 . Taylor and Randeria 72 established sum rules at unitarity for the dynam- 
ical counterpart rj(ui) and £2(0*), which are less general than those in Eqs. (18) and (39). Zwerger and colleagues 73 addressed 
the behavior of the normal unitary gas predicting a minimum in 77/s near the superfluid transition. 

Bruun and Smith 47 ' 74 were, perhaps, the first to emphasize the importance of the (static) shear viscosity in cold Fermi gases 
using both a high temperature fermionic Boltzmann approach and a Kubo-stress tensor based scheme (above T c ), within BCS- 
BEC theory. They importantly recognized 47 that the introduction of a pseudogap would lower the normal state 77. However, the 
diagram set which was used was "not conserving" 47 . 

Rupak and Schafer 4 introduced an alternative (bosonic Boltzmann) transport theory in which the low T «T C shear viscosity 
is dominated by the Goldstone bosons or phonons. They predicted that 77/s, increases as temperature decreases. In their approach 
phonons dominate the low T transport, presumably leading to the same upturn in 77, as is seen in helium-4, and shown in Figure 
1, on the right. Establishing diagrammatic consistency is a central theme of this paper and within a BCS-like theory the shear 
viscosity, as a transverse response, does not directly couple 16 to the Nambu-Goldstone bosons. This is to be distinguished from 
the longitudinal current-current correlations which are the basis for the bulk viscosity and necessarily depend on these collective 
modes. 



C. Transport in High Tc superconductors 

Understanding the temperature dependence of the (inverse) conductivity or dc resistivity, particularly near optimal doping, 
was one of the first puzzles posed by the high temperature superconductors. It should, nevertheless be noted that a number of 
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ideas and interpretations about transport were established well before there was a full understanding of how widespread was the 
pseudogap phase (which extends over most of the phase diagram, including optimal doping). For this reason, one can argue that 
these earlier theories must be reexamined. 

A number of different models from spin-charge separation 66 to marginal Fermi liquid phenomenology 75 were invoked to 
explain (i) the unusual normal state temperature power laws observed, e.g., the linear resistivity. A second notable observation 
(ii) was that the effective carrier number in transport appeared to scale with the number of (extra) doped holes 66 , called x, even 
though the volume of the Fermi surface scaled as expected with 1 + x. The first of these experimental observations gave rise 
to the concept of a distinct "strange metal" phase which was thought to separate the overdoped and underdoped regimes. The 
second of these gave rise to the concept of "spin-charge" separation. Despite the widespread agreement that the normal state of 
the cuprates was a non-Fermi liquid phase, there is today a strong belief that the superconducting phase is Fermi liquid based 66 . 

It must be noted that the BCS-BEC crossover scenario addresses these issues from a very different perspective. It presumes that 
there is a smooth evolution, as seen experimentally 1 1 from over to underdoped behavior rather than a distinct strange metal phase 
near optimal doping. We have argued here that the suppression in the number of carriers may be associated with the fact that the 
magnitude of the excitation gap grows with underdoping. [In the ac conductivity, this leads to a very small weight for the ui w 
Drude peak, which requires by the sum rule, an additional mid-infrared contribution, now widely observed experimentally 21 .] 
As noted above, there is no spin-charge separation in the normal phase. Importantly, the superconducting phase has bosonic 
excitations representing non-condensed pairs (also present above T c ) and therefore this phase is not Fermi liquid based. Finally, 
the dissipation mechanism here associated with the BCS-BEC crossover scenario arises from the inter-conversion of fermions 
and pairs, which is distinctly tied to the pseudogap. This is to be distinguished from most transport theories 76 77 which focus on 
impurity effects or the dynamics of the pairing boson. 

There is a shortcoming in our cuprate transport calculations because we have ignored impurity effects in transport altogether. 
Away from very low T, it is generally accepted 21 that they are not particularly important. Nevertheless, they may lead to 
incomplete condensation 30 ' 77 in the ground state. Indeed, it is difficult to see how to reconcile claims 21 that the T = superfluid 
density scales as oc x, with the transverse f-sum rule (Eq. 39)), without invoking incomplete condensation in the ground state. 
A deeper understanding of a possibly more inhomogeneous treatment of impurities is needed. Fortunately, this is not an issue in 
the cold gases. 

With a growing appreciation for the nature of the pseudogap (and related "Fermi arc" effects), experimentalists have provided 
some support to our findings. From Ref.78 it is said that "This indicates that the functional form of the dc resistivity of cuprates 
Pdc(T) is governed not only by the relaxation processes but also by temperature-dependent numbers of carriers". Moreover, in 
Ref.79 it is stated that "One may notice that a natural extension of the present argument would be that the T-linear resistivity 
usually observed near optimum doping may not necessarily by a sign of a T-linear relaxation rate, because n e ff may be changing 
with T". An interesting corollary to a carrier number which necessary increases with T, is that the inverse lifetime 7 should 
contain higher powers than linear (we use the most conventional, Fermi liquid dependence 7 oc T 2 ) to arrive at metallic behavior 
for the conductivity. 

VII. CONCLUSIONS 

In summary, this paper has addressed the role of the pseudogap in the oj — > conductivity and in the shear viscosity both above 
and below T c . We have emphasized the analogy between "bad metals" 2 and "perfect fluids" 1 seen in high T c superconductors 
and the atomic Fermi gases. Both of these phenomena, we argue, may arise from pseudogap effects. 

Our approach builds on a consistent gauge invariant treatment of transport (which has not been addressed previously), in 
which the transverse f- sum rule demonstrably holds. In this paper we have demonstrated success in simultaneously addressing 
experiments in cold gases and high T c cuprates within the same transport formalism. As a summary, Figure 5 gives a reasonable 
understanding of the experiments shown in Figures 3 and 2. Equally important, Figure 6 shows semi-quantitative agreement 
with shear viscosity data from below to above T c . We see no sign of the upturn which others have predicted. 

It was our intention in this overview to introduce some of the key challenges in understanding transport in the high T c 
superconductors to the wider readership interested in perfect fluidity. It is hoped that an appreciation of this broader context may 
lead to new breakthroughs in understanding superficially distinct, but quite possibly connected, physical systems. 
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Appendix A: Rewriting the number equation 

It is useful to fill in a few steps in the line of argumentation from the text. We wish to show here for the simpler case of a 
Fermi gas, how the number equation may be rewritten using a Ward identity: 

P(Q) = E( fc + |) 2 [ G kGk + q + t P9 (P)Go,p-K- Q Go,P-KG K+Q G K } (Al) 

K P 



And the number equation 



n 



K K 

= 2Y,k a G 2 (K)dG K 1 /dk a 

K 

= 2 ]T k a G 2 K [dG-^/dK - dx(K)/dk 



K 



= -2j2k a Gl[^-J2t(Q)dG^ Q _ K /dk a ] 



m 

K Q 



n 



K Q m 

= -\T, k ^[ G K+Y, t ^Q)Gl,Q-KG 2 K ] (A2) 

O Tib „ 

K Q 



where the last step specializes the result to T > T c . The last equation is the central result we wanted to prove: from which we 
have 

2 

716 

K(0) = + P(0) = 0, above T c (A3) 

m 

Appendix B: NSR Theory and the Normal state Structure Factor 

In a recent Physical Review Letter 14 we have shown how two photon Bragg scattering can be used to establish in situ the 
presence of coherent order in a superfiuid, at any temperature, wavevector and frequency. For the most part experiments on 
unitary gases have relied on sweeps to the BEC, to find evidence for condensation. Our analysis is based on the definitions and 
sum rule in Eq. (24). It, thus, depends on imposing the current conservation laws, which have been extensively studied and 
verified 14 . 

As an alternate example, here we consider Nozieres Schmitt-Rink (NSR) theory 35 in the normal state. Figure 7 indicates the 
characteristic class of diagrams. These diagrams enter into the density-density and spin-spin correlation functions. The calcu- 
lations for the spin response build on an earlier publication 17 . We define the pair susceptibility Xo(Q) = Cot,pGo4,.Q-p 
where Go<j is the bare Green's function. The pair propagator is given by t Q (Q) = U/(l + U\o(Q)), where U is the two- 
body interaction. Here we present a more systematic discussion of the spin response than what was shown in Sec. III. The 
interaction term in the Hamiltonian is given by if / ~ / d 3 rJ^A s ^ where A$ = (B z , m) is the "effective" 4-vector field and 
= (n s , J s ). Here B z is the z component of the magnetic field, m is the magnetizing field, n s is the z component of spin 
and J s is the magnetization current. The spin rotational symmetry around the z axis leads to the conservation law of spin: 

Jf = 0. 

From linear response theory the spin response kernel can be written as 



QlAQ) =J2^AP,P + Q)Goa,p+QAL(P + Q,P)G 0a ,P + ^9^(1-9^), (Bl) 



where A^ CT (F, P + Q) = (^(p + %),S a ) is the bare vertex function of the spin-external field interaction, and is the full 
vertex function. There is an implicit summation over the indices a. Note, importantly, that the vertex function has different 
signs for different spin indices, which is to be contrasted with the charge or equivalently density response functions. In order to 
satisfy local conservation laws, the vertex must satisfy a Ward identity: Q ■ J2p ^%(Pi P + Q) = Sa{pQ a ,p ~ ^OctP+q)- As 
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a result, in Nozieres Schmitt-Rink theory there are three contributions = + <JA^ MT + <5A^ AL , where the subscript MT is 
associated with the contribution from the Maki-Thompson (MT) like diagrams and AL the Aslamazov-Larkin (AL) diagrams. 
An important result (for singlet pairing) is that the contribution from the AL diagrams automatically vanishes. Thus, for the spin 

structure factor Ss = Sso + S MT s with S MT s = —S MT c. 

pg pg pg 

By contrast, in the normal state, the dynamical structure factor in NSR theory, for the particle density 17 can be written as the 
sum Sc — Sco + S MT c + 2Sal- Here one should note that Sso — Sco- The remaining terms denote the corrections from the 
MT, and AL diagrams. Note that the spin and charge contributions from the MT term enter with opposite signs. Moreover, the 

Ward identity implies a cancellation 6 S MT c + 25 , A l = —S MT c. Thus 

pg pg 

Ss = Sso + >Sm T s = Sc = Sco - S MT c, (B2) 

which proves the desired result: in the normal state the spin and charge degrees of freedom are indistinguishable. Note that 
the theoretical proof of this result depends on using a consistent theory of BCS-BEC crossover with full gauge invariance. The 
same diagrams that are used in the above proof are needed to show that there is no Meissner effect in the normal state. In this 
way we have demonstrated 

5_(cj,q)=0 above the transition. (B3) 

We are now led to an important observation: the quantity S- (u>, q) for all q, uj, 1 /kpa can be used as an indication of in- situ 
superfluid order, without requiring sweeps to the BEC. We show a plot of this behavior in the right side of Figure 7 where it can 
be seen that the difference structure factor vanishes in the normal state. The behavior below T c is shown in the figure for the 
case of BCS-Leggett theory 17 . 
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